Average persistence in random walks 
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Abstract. — We study the first passage time properties of an integrated Brownian curve 
both in homogeneous and disordered environments. In a disordered medium we relate the 
scahng properties of this center of mass persistence of a random walker to the average 
persistence, the latter being the probability Ppr(i) that the expectation value of the 

walker's position after time t has not returned to the initial value. The average persistence 
is then connected to the statistics of extreme events of homogeneous random walks which 
can be computed exactly for moderate system sizes. As a result we obtain a logarith- 
mic dependence Ppr{t) ~ ln(t)^^ with a new exponent 6 = 0.191 ± 0.002. We note on a 
complete correspondence between the average persistence of random walks and the mag- 
netization autocorrelation function of the transverse- field Ising chain, in the homogeneous 
and disordered case. 



First passage time or persistence problems have a long history in the physical literature 
[jl|. Recently they gained a lot of interest, since the persistence exponents, describing the 
asymptotic behavior of first passage time probabilities, are shown to be independent dynami- 
cal critical exponents, which have been calculated for various models exactly]^, Not much 
is known about analogous quantities in systems with quenched disorder, for instance the ran- 
dom walk (or diffusion) in a disordered environment @] , which is in the one-dimensional case 
the Sinai- model 0. For this model the first passage / persistence exponent for a single 
walker has been determined by us in a previous work ^j. In this letter we introduce and 
study the concept of average persistence of random walks both in homogeneous and random 
environments. 

We consider a random walk with nearest neighbor hopping in one dimension defined by 
the Master equation 

Piit = 0) = 5i^i , ^Pi(^) = -{wi^i^i + Wi^i+i)pi{t) + Wi_i,iPi_i(t) + Wi+isPi+iit) (1) 

describing the time evolution of the probability Pi{t) for the walker to be at site i after 
time t when having been initially at site i = 1. The homogeneous random walk is defined 
via uniform transition rates = Wj+i^j = 1/2 and the random walk in a disordered 

environment is modeled by choosing the transition rates to be quenched random variables 
that obey a particular distribution, e.g. the uniform distribution vTuni given by Wi^i±i = 1 for 
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Figure 1: The average persistence [Ppr(L,t)]av 
to the power —1/29 a versus Int, which should 
yield a straight line in the limit L — > cxd, for 
different finite size scaling length L. The data 
are for uniform distribution and averaged over ^ 
50.000 samples. We used an estimate 6a = 0.19 
to obtain an optimal data collapse in the scaling J 
region. The straight line indicates the asymp- o! 
totic result for the infinite system L — > oo. 
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< Wi^i±i < 1 and otherwise, or the binary distribution 

1 1 

T^hin{Wi^i-l) = -5{Wi^i-l- X) + -5{Wi^i^l- \'^) , = 1 Vi (2) 

with A some arbitrary parameter. For the disordered case physical observables have to be 
averaged over this distribution tt which is denoted by square brackets [. . .]av Note we consider 
the general case of asymmetric hopping rates ti^j^j+i ^ t^i+i^j and we do not confine ourselves 
to the so called random force model with correlated transition probabilities parameterized 
as Wi^i±i = Aexp{±(f)i) with random, uncorrelated potentials (pi on each site. 

In order to define single walker persistence probabilities we put an adsorbing boundary 
at site i = 0, which means that we set Wi^o = and also introduce a finite size length 
scale L into the system by putting another adsorbing boundary at z = L + 1, i.e. setting 
wl+i,l = 0. Now we define the length scale dependent single walker persistence Ppr(-L,t) to 
be the probability that a walker does not cross its starting point (i.e. does not get trapped at 
site i = 0) within the time interval t. The following arguments[^] will then lead us to a scaling 
form for Ppj.{L,t): a) the typical time the walker needs to reach the site L + 1 scales like 
ttyp ~ in the case with symmetric transition rates and like In t^yp ~ VL in the asymmetric 
case, the Sinai- model, b) the asymptotic limit Psurv(-^) = hnij^oo Pprit, L) = \imt^ooPL+i(t), 
which is what we call a survival probability, behaves like Psurv(-^) ~ in the symmetric 
case and like Psurv(-^) ~ L"^" with 6 = 1/2 in the Sinai-model. Hence we expect 

/ L~^' ■ps{L'^lt) (symmetric case) 

[^prli^, rjjav ~ j j^-e^ . Pa{^/ In t) (asymmetric case) ^ ' 

where the persistence exponents are = 1 for the symmetric case and 9a = 1/2 for the 
asymmetric case, and the scaling functions behave like ps{x) x^"^"^ and pa{x) x^^" for 
X — *■ CX3 and ps/a{x) — >■ const, for a; — *• 0. In the infinite system size limit one thus has for 
persistence probability Ppr(i) in the asymmetric case Ppr(t) = liniL^oo Ppr{L, t) oc ln(t)~^^'', a 
logarithmically slow decay reminiscent of the critical dynamics of the surface magnetization 
of the random transverse Ising chain[0] (RTIC). This is not incidental: the equivalence of 



the surface magnetization in the latter quantum spin chain and the survival probabilities 
of random walks has been formulated the first time in ^ and further analogies between 
anomalous diffusion and the RTIC have been uncovered recently [0, |]. 

It is known that diffusion in the Sinai model is different from normal diffusion in many 
respects[ffl. Consider for instance one particular disorder realization. Then an initially 
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narrow probability distribution of a walker peaked around, say, x(t = 0) = does not 
broaden with time, only its expectation value {x{t)) = ■ Pi{t) diffuses logarithmically 
slowly away from its starting point. Therefore, in this situation in addition to the single 
walker persistence one should also consider the persistence properties of the average position 
of the walker the average persistence Ppr(t): This is the probability that up to a 

specified time t the average position of the walk in a particular environment has always been 
on one side of the starting point, i.e. \/0 < t' < t : {x{t')) > 0. 

To study the average persistence we consider again a finite size situation where we put 
an adsorbing boundary at i = —L and one at i = L (note that now a single walker is not 
adsorbed when crossing the starting point but only when he leaves the finite strip of width 
2L centered around i = 0). Then the average persistence Ppr(L, t) obeys the same scaling 
form as in the second line of @, but the single walker persistence exponent 9a replaced by 
the average persistence exponent 6 a- 

[Pp,(L,t)]av~L-^"-p„(v^/lnt). (4) 

Moreover, the persistence probability for the infinite system decays again logarithmically 
slow: Ppr{t) oc ln(t)~^^''. Recently it has been conjectured for the random force modelf^ 
that the exponent 6a is related to the golden mean via ^„ = (3 - v/5)/4 ^ 0.191. In fig. 1 
we show numerical data that have been obtained by a numerical calculation of PpT-{L,t) via 
diagonalization of the linear operator on the r.h.s. of eq.(P, which indicate that this might 
also hold for the general asymmetric case. 

In what follows we will demonstrate how to obtain a precise estimate of the average 
persistence exponent 9a- Since in the limit t ^ oo the average persistence probability is given 
by Psurv ~ L'^"-, the computation of the survival probability Psurv will lead us to the desired 
result. We will now relate this survival probability to a problem in the statistics of extreme 
events of homogeneous random walks. To this end we consider the binary distribution vr in 
d^) in the limit A — *• 0, which means that one can discriminate between "forward" bonds 
+ 1), which are those for which = A — > 0, i.e. those where almost only (with 

probability 1/(1 + A) —>■ 1 jumps from i to i + 1 occur, and "backward" bonds, which have 
Wi^i-i = A^^ oo, implying almost always jumps from i to i — 1 (again with probability 
A~^/(l + A"^^) 1). This implies that when sketching the transition rates as configuration 
as is done in fig. 2 the disorder configuration can be done visualized as a random landscape 
with hills and valleys. Thus a walker starting at i = is on average first driven to the 
first minimum Mi, and considering the finite size situation where L = Li this disorder 
configuration counts as a surviving configuration, since the walker in average spends most of 
its time at Mi. Increasing the length scale L further beyond L2 will then drive the average 
position of the walker over the barrier B into the valley M3 at negative coordinates, which 
means that this configuration is dead, i.e not surviving. 

Thus we conclude that the computation of Psurv amounts to counting all surviving config- 
urations of a random landscape, i.e. a homogeneous random walk, in the manner described 
above and the ratio of surviving configurations is just Psurv In other words we consider 
random walks of length 2L, corresponding to realizations of the transition probabilities ac- 
cording to (1^) in a strip of width 2L centered around the starting point, and say that it 
is in average surviving if certain conditions are fulfilled. These conditions are checked by 
inspecting the random landscape generated by the disorder configuration (i.e the transition 
rates): one scans the landscape in both directions from the starting point and denotes with 
h{i) {i = —L, —L + 1, . . . , —1, 0, 1, 2, . . . , L) the position of the walker (or height of the land- 
scape) at step (or site) i. We define the extreme events to the right and to the left of z = 
by Xniax/min(0 = max/mm{h{j)\0 < j < i} and ?/max/mm(0 = max/min{/i(j)| - i < j <0}, 



Figure 2: Sketch of the random landscape that 
corresponds to the binary distribution (^). In 
the hmit A — > transitions corresponding to 
broken arrows will rarely occur whereas tran- 
sitions corresponding to full arrows will occur 
with probability close to one. For further de- 
tails describing the effective dynamics see text. 
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respectively, and check iteratively (from i = 1 to i = L) whether 

(i) and Xmax(«) > l/max(0 • (5) 

If this happens for some site as it does for i = 8 in fig. 2, it means that there is a lower 
minimum on the left side of the starting point (xmin < 1/min) and that the walker can go there 
since the barrier in between is low enough (xmax < 2/max)- This implies that this configuration 
is dead. 

In the inset of fig. 3 we show the results for a numerical estimate of the survival probability 
Psmv{L) by inspecting 10^ random walk configuration for different system sizes — the data 
fit well to a power law with the exponent 6a = 0.19. Next we implemented a recursive routine 
that computes the number of surviving configuration according to the above extreme events 
criterion (|^) exactly. Here we took special care to the degenerate minima (like Mi and M2 
in fig. 2), in which case we assumed the arithmetic mean of the location of the two to be the 
effective position (Li in fig. 2). From these exact data for system sizes up to L = 14 we can 
extract via Oa{La, Lb) = — \n{Psnrv{La) / Psmv{Li,)) / \ia{La/ Lb) an effective finite size exponent 
that approaches the exact value. In fig. 3 we show these exact data for La = L = Lb + 1 
and La = L = Lb + 2. Finally, from a standard series extrapolation procedure0 applied to 
these data we obtain our best estimate which is 

= 0.191 ±0.002 , (6) 

in very good agreement with the conjectured value = (3 — v^)/4 = 0.19098. . . for the 
random force model 0. 

Although the concept of an average persistence seems to be based upon the special 
feature of non-dispersive diffusion in the Sinai model, which we described above, it turns 
out that there is an analogy to it in normal diffusion, too. To show this, we first introduce 
the integrated position of the walker at step t a.s It = Yli=i V, where v is its position at 
time-step r. Then, we define the center of mass persistence through the survival condition 
for the integrated position as If > for < t' < t. For the disordered case we define a 
center of mass position via 

C{z) =±j- M{h{j)) . (7) 
j=-i 

where the weights M{h{j)), which are proportional to the average time the walker spent at 
a given site, depend on the random configurations. In a strip (— L < i < L) with adsorbing 
boundaries the survival condition of the average position in the large t-limit is C{i) > for 
i = 1, . . . ,L, which is equivalent to the survival condition in (|^). For the extreme binary 
distribution, we consider here, the weight function is singular, such that M{h')/M{h) 
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Figure 3: The effective finite size exponent 9a 
from exact values for Pp^{L) via complete enu- 
meration of all surviving configuration. The 
asymptotic value 6a = 0.191 obtained by se- 
ries extrapolation methods is also shown. The 
insert shows the result for PpT-{L) via stochas- 
tic enumeration for larger system sizes, giving 
da = 0.19 ±0.01. 
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for h' < h, which is a consequence of the hmit A — *• for the distribution (|^). Since the form 
of the distribution of the random transition rates is generally irrelevant one expects that the 
scaling properties of the average persistence and those of the center of mass persistence in 
d^) are equivalent for other type of distributions as well. 

For a homogeneous walk one can also study the center of mass persistence or the survival 
probability of the If integrated position. This type of problem has already been considered by 
Sinaifl^ for a discrete model and later by Burkhardt[0 for the continuum model. According 



to these exact results in a homogeneous infinite system the long-time behavior of the average 
persistence is given by hmL~^oo Ppr iL,t) ~ t~^^^. We use this result to write for finite 
systems the scaling conjecture: 

PTiL,t) = L-~'^-p,,Mt) (8) 

with 9 horn = 1/4 and the scaling function behaves like Phomi^) ~ x^'*"™ for a; — *• cxo and 
Phomi^) ~^ const ioT X —> 0. The scaling combination, x = L/t in follows from the scaling 
properties of the Brownian motion. In t steps the walker typically crosses its starting position 
~ t^/^ times and its displacement is typically of ~ t^/^, thus the integrated displacement 
is ~ t, from which the scaling relation L ^ t follows. It is interesting to note that the 
integrated Brownian curve is an isotropic critical object, with a correlation length critical 
exponent z/ = 1, as for the two-dimensional Ising model. 

By an exact enumeration up to L = 30 we have checked the relation in (H) in the limit 
t oo and got an estimate p^omi.^ = 0) = 0.6183(3) for the scaling function. It is amusing 
to note that the above value agrees, within the accuracy of the estimate, with the inverse of 
the golden mean ratio 1/r = 2/(^5 + 1) = 0.6180. 

Concluding, we introduced the concept of average persistence in random walks and stud- 
ied its scaling properties in homogeneous and in disordered environments. The finite-size 
scaling behavior of Ppr{L,t) as given in and (H) involve new exponents: 6 in (j^) and 
dhom = 1/4 for the disordered and homogeneous problems, respectively. It is interesting to 
note on an analogy between the finite-size scaling form of the average persistence of random 
walks and that of the magnetization autocorrelation function G{L,t) = [(o'£/2(0)cr£/2(^)]af 
of the transverse-field Ising spin chain of length L. According to exact [p!^ and conjectured 
plus numerical ||13|| results we have for the finite-size scaling behavior ||1 4] | 



G(Lt)r^l ^ 3hom{L/t) (homogeneous case) 

' ■ ^j.a„(i(\/L/ Int) (random case) 




with the magnetization scahng dimensions Xhom = 1/8 and Xrand = (3 — V^)/2 ~ 0.191||T3|, 
for the homogeneous and random transverse-field Ising models, respectively. Since 2xhom = 
9hom and Xrand = Oa we havc a complete correspondence between the scaling behavior of the 
average persistence of random walks and the magnetization autocorrelation function of the 
transverse-field Ising chain, both in the homogeneous and in the disordered problems. These 
relations then complete the previously observed correspondences between single walker 
persistence and the surface autocorrelation function of the transverse-field Ising spin chain. 
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